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Abstract. The infrared (IR) emission plays a crucial role in understanding the star formation in galaxie s hidden by 
dust. W e first examined four esti mators of the IR lumi nosity of galaxies, Lfir (IHelou et al Lttr IIDale et alJ 

l2001al) . revised version of Ltir ijDale fc Heloull2002h (we denote Li-ma), and Ltr (|Sanders fc MirabeJIlQQdT bv 
using the observed SEDs of well-known galaxies. We found that Lir provides excellent estimates of the total IR 
luminosity for a variety of galaxy SEDs. The performance of L tir2 was also found t o be very good. Using Lir, 
we then statistically analyzed the IRAS PSCz galaxy sample l|Saunders et al.l|200(tl and found useful formulae 
relating the MIR monochromatic luminosities [L(12 pm) and L(25 pm)] and Lir. For this purpose we constructed 
a subsample of 1420 galaxies with all four IRAS band (12, 25, 60, and 100 pm) flux densities. We found linear 
relations between Lir and MIR luminosities, L(12 pm) and L(25 pm). The prediction error with a 95 % confidence 
level is a factor of 4-5. Hence, these formulae are useful for the estimation of the total IR luminosity only from 
12 pm or 25 /im observations. We further tried to make an 'interpolation' formula for galaxies at < z < 1. For 
this purpose we construct the formula of the relation between 15-pm luminosity and the total IR luminosity. We 
conclude that the 15-pm formula can be used as an estimator of the total IR luminosity from 24 pm observation 
of galaxies at z ~ 0.6. 
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1. Introduction 

Star formation activity is one of the fundamental prop- 
erties useful to explore the evolution of galaxies in the 
universe. Generally, the star formation rate is measured 
by the emission from young stars, i.e., ultraviolet (UV) 
and related nebular line emissions. However, a significant 
fraction of UV photons are absorbed and re-emitted by 
dust mainly in the infrared (IR), hence the IR emission 
plays a crucial role for an unders tanding of the obscured 
star formation in ga laxies (e.g., iBuat et aTlll999l 120021: 
iHirashita et~alll2003fl . 

Further, clarifying the correlation between flux den- 
sities at various IR bands is an important task to under- 
stand the origin, release and transfer of energy in galaxies. 
Such studies play a crucial role in constr ucting and ver- 
ifying IR galaxy evolution mo d els (e.g.. iGran ato et all 
200ri iFranceschini et all l200lt iTakeuchi et all l2001alrl 



Takaei et alJl2003|) . 



Send offprint requests to: T. T. Takeuchi 
e-mail: tsutomu.takeuchi@oamp.fr. 

* Postdoctoral Fellow of the Japan Society for the Promotion 
of Science for Research Abroad. 



Based on their 12-pm sample of galaxies, 
ISpinoglio etafl l|l995|) made a pioneering study to 
examine various correlations between flux densities 
from near-IR (NIR) to far-IR (FIR), and presented 
useful diagnostics for Seyferts and normal galaxies on 
color-color diagrams. They also found that the 12-pm lu- 
minosity correlates well with the bolometric (0.4-300 pm) 
luminosity. 

Now that data obtained by Spitzer have started to be- 
come available, we are better able to explore the IR prop- 
erties of galaxies at high redshift. 1 The 24- um. band o f 
Spitzer MIPS is very sensitive (e.g.. lPapovich et afll2004h . 
and will be used extensively for the studies of high-z galax- 
ies. Hence, from a practical point of view, it is worthwhile 
to find a good estimation method of the total IR luminos- 
ity of galaxies from the mid-IR (MIR) luminosity. This 
will also be useful for forthcoming IR space missions, e.g., 
ASTRO-F. 2 

In this work, we present the estimation formulae for 
the FIR luminosity from the MIR. We focus on the rela- 
tion between MIR and total IR luminosities, in contrast to 



1 URL: http://www.spitzer.caltech.edu/ 
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ISpinoglio et alJ (^995) , who used the bolometric luminos- 
ity integrated from the optical to the IR. For this purpose, 
we have to rely on some conventional formulae to estimate 
the total IR luminosity, since direct measurement of the 
total IR luminosity is possible only for a limited number 
of galaxies. First we examine the performance of four for- 
mulae in use, using galaxies with well-measured spectral 
energy distributions (SEDs). This sample consists of 17 
galaxies ranging from dwarfs to ultraluminous, and from 
cool (submillimetrc bright) to hot (MIR bright) ones. 

We then perform a correlation analysis for the galaxy 
sample extracted from IRAS PSCz, and obtain a statisti- 
cal formula for the estimation of the total IR luminosity 
from MIR luminosities. This statistical sample is selected 
by the criterion that the galaxy has all four IRAS flux 
density values. By combining the formula and ISOCAM 
15-/mi data, we then give an interpolation formula of the 
FIR luminosity for galaxies at z ~ 0.6 observed in the 
Spitzer MIPS 24-/xm band. 

The paper is organized as follows: we examine the four 
estimators of the total IR luminosity in Sect.|2 We present 
our statistical analysis based on IRAS PSCz galaxies in 
Sect. El A reexamination of the estimator Lir and ap- 
plication to galaxies at z ~ 0.6 are given in Sect. 0] 
Sect. [S] is devoted to our conclusions. The SEDs of ob- 
served galaxies used in Sect. [5] are shown in Appendix A. 
Mathematical details of the regression analysis are pre- 
sented in Appendix B. 

We denote the flux densities at a wavelength A by a 
symbol S"a, but the unit is [Jy]. Throughout this work, we 
assume a flat lambda-dominated low-density universe with 
cosmological parameter set (h, Q,q,Xo) = (0.7,0.3,0.7), 
where h = Ly 100 [kms -1 Mpc _1 ]. 

2. Performance of the estimators for the total IR 
luminosity 

Since direct measurement of the total IR luminosity is 
only available for a limited number of galaxies, we have 
to use a formula to estimate the total IR luminosity from 
discrete photometric data, mainly in the IRAS bands. In 
this section, we examine the performance of four formulae 
in use. 

2.1. Estimators 

First, we define £„ as the luminosity per unit frequency 
at a frequency v = c/A (c: the speed of light). The unit of 
C v is [ergs _1 Hz _1 ] throughout this work. 

We examine the following four total IR luminosity es- 
timators. 

1. The classical FIR luminosity between A = 42-122 /im 
(jHelou et al.lll988j) . defined as 

L Fm = 3.29 x 10~ 22 

x (2.58/^(60/^1) + C(100/im)) [L ] . (1) 



Table 1. Well-known galaxy sample. 



Name 


References' 1 


Normal galaxies (10 9 Lq 


< L m talb < lO n L ) 


M 63 


1 


M 66 


1 


M 82 


1 


M 83 


1 


NGC 891 


1 


NGC 3079 


1 


NGC 4418 


1 


NGC 7714 


1 


IR luminous galaxies 


(L m tal > lO 11 !^) 


NGC 2623 


1,2 


NGC 7679 


1,2 


UGC 2982 


1,2 


UGC 8387 


1,2 


Arp 220 


1,3 


IRAS F10214+4724 


1,4 


Dwarf galaxies (L 


m tal < W' J L e ) 


NGC 1569 


1,5 


II Zw 40 c 


1,6 


SBS 0335-052 c 


1,7,8 



a References: 1. NED , 2. iDunne fc Ealed (l200lt) . 3. 
Downes et"aT1 Jl993f) . 4. iDownes et alJ lll992l). 5. 
Galliano et all ll2003h. 6.lMaddenl (I2000T) . and 7. bale et all 
feOOlbl) . 8. lHouck et all (|2004) . 

b Total IR luminosity Lir 3,1 is calculated by integrating over 

the wavelength range of A = 8-1000 /tm. 
c The longest wavelength flux den sities are calcula t ed by an 

extrapolation using the model of lTakeuchi et al.l (l2003al) . 



2. The 'total' IR luminosity, TIR (A = 3 -1100 /im), ad- 
vocated bv lDale et all l|200l4 

L TIR EE L F IR X IQ^o+^x+a^+asx'+a^ 4 ^ ^ 

where x = log(S , 60 /S'ioo), and (a , a x , a 2 , a 3 , a 4 ) = 
(0.2738, -0.0282, 0.7281, 0.6208, 0.9118). 

3. The upd ated version of t he TI R (A = 3-1100 /mi) pre- 
sented bv lDale fc Heloul l)2002(l (here we denote Ltirs) 

L T ir2 = 2.403 vC v {2h^m) - 0.2454i>£„(60/im) 

+1.6381^(100/im) [Lq]. (3) 

This is better calibrated for submillimeter wavelengths 
than Ltir- 

4. The luminosity between A = 8-1000 /mi presented by 
ISanders fc Mirabel ([i.996). In this work we refer to 
their IR luminosity estimator as Lir. 

Lir = 4.93 x 10~ 22 [13.48£„(12jum) + 5.16A,(25/rai) 
+2.58£„(60/mi) + C u (l00fim)} [Lq] . (4) 

2.2. Examination of the IR luminosity estimators by 
known galaxies 

Though these IR luminosity estimators are popular in 
related fields, direct comparison between the measured 
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Fig. 1. Comparison between the true total infrared 
luminosity and the estimated luminosity of galaxies. 
Open triangles, open diamonds, open squares and filled 
squares represent the ratios Lpm/L^ 3,1 , LTiR./£iR tal , 
£ T iR 2 /£& tal , and L m /L\^\ respectively. 



IR luminosity and the estimated value has rarely been 
done to date. We examine the performance of the above 
estimators using the SEDs of observed galaxies. We 
compiled 17 galaxies with well-measured flux densities, 
with a total IR luminosity range of 10 8 L© ^ £iR tal ^ 
10 13 L© (see Table QJ. Among the dwarf galaxy sam- 



ple (L\$ ai < 10 9 £ Q ), the longest wavelength data (i.e., 
FIR and submm) are not available for II Zw 40 and 
SBS 0335—052. We calculated the flux densities by ex- 
trapola ting their SEDs using the model o f lTakeuchi et alJ 
l|2003a|) (see also iTakeuchi fc Ishiil EoO^ . The compiled 
SEDs are presented in Appendix A. 

We calculated £jR tal by integrating the observed data 
directly within a wavelength range of A = 8-1000 /im 
by interpolation and extrapolation. Figure ^ shows the 
comparison between L\^ al and estimated luminosities of 
galaxies. As expected, the classical Lfir gives system- 
atically lower luminosities than the true ones, because 
it represents the luminosity at 42-122 /im, and there- 
fore the MIR and submm radiations are not included. 
Especially, two galaxies with hot dust (SBS 0335—052 and 
IRAS F10214+4724) significantly deviate downward from 
t h e diagonal line. 

bale et ail l|2001a|) considered the correction factor for 
the contribution outside the range of Lfir as a function 
of the ratio S^o/Sioc)- We see that the estimation is clearly 
improved, but the IR luminosities of the two extreme ob- 
jects are still underestimated. This is because their Ltir 
has been designed for normal galaxies, and not for such 
extreme objects. 



In contrast to the above two estimators, £tir2 and 
Lin give much better estimates for all the galaxies in 
Tabled They work not only for the objects with very hot 
dust emission like IRAS F10214+4724 and SBS 0335-052, 
but also for a heavily extinguished galaxy like Arp 220. 
For SBS 0335—052, Lift gives a better result. This is an 
expected result because Ltir2 uses three (25, 60, and 
100 /im), and Lift uses all four IRAS flux densities. In 
general, Ltir2 gives slightly larger values than Lift does, 
probably because the considered wavelength range for the 
former (A = 3-1100 /im) is wider than that for the latter 
(A = 8-1000 urn). 

Thus, Lift is the best estimator of the total IR lumi- 
nosity. As long as we have the four IRAS flux densities, we 
can obtain a precise estimate for the total IR luminosity. 
When data in three (25, 60 and 100 /im) or two (60 and 
100 fim) bands are available, Ltir2 and Ltir give reason- 
able values except for galaxies with extremely hot dust. 
Ltir2 works almost as accurately as Lift. In the following 
discussions we regard Lift as the correct estimate of Lj^ tal 
and use Lift as Lj R tal itself. 

3. Statistical analysis of the IRAS sample 

Our next step is to find a conventional formula to estimate 
Lift only from a single MIR band. For this purpose, we 
make a regression analysis for Lift and MIR luminosities 
in the IRAS bands. Here we define the luminosity at a 
wavelength A, L(X), as 



L(A) = XC X = vL v 



(5) 



and we discuss L(12/im) and i(25/im). Mathematical 
details of the regression analysis can be found in 
Appendix B. 

3.1. IRAS Sample 

We selected a sample from IRAS PSCz (hereafter PSCz, 
ISaunders etHI kdOft . The PSCz is a complete, flux- 
limited all-sky rcdshift survey catalog of IRAS galaxies 
with a detection limit of Sqo > 0.6 Jy. It contains 15411 
IRAS galaxies with redshifts. Out of the whole sample, 
we selected galaxies with good quality flux densities for all 
four IRAS bands (12, 25, 60, and 100 /im) for this analysis, 
because Lift requires all four flux densities. We performed 
this procedure as follows: 1. We examined the flux origin 
and quality flags given in PSCz for the point source flux 
density, and omitted galaxies with upper limits (denoted 
as 1 in pscz.dat), 2. We extracted the coadded or extended 
addscan flux densities. We adopted this selection because 
we found that the addscan/coadded fluxes with quality 
flag 1 include unrealistic values close to the upper limits 
in point source flux densities. 

There is a caveat that the selection by using all the 
four IRAS bands would introduce a subtle sample bias in 
the analysis. In order to see if the bias is serious, we also 
made a subsample by omitting the galaxies with flag 1 only 
at 100 fim (3260 galaxies included). This subsample for 
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Fig. 2. The relation between 12 /mi luminosity L(12 /im) 
and Lir, the estim ate of the total infra r ed (IR ) luminosity 
L totai p r0 p 0sec i bv lSanders k Mirabell l)l996h . 

comparison gave essentially the same result as the above 
sample (the difference was less than ~ 1 %). This means 
that the selection in the MIR affects the result only very 
slightly, and the sample properties are c ontrolled by the 
FIR. I t is a clear contrast to the sample of lSpinoglio et all 
l|l995h which was 12-/mi selected: the present sample con- 
sists of more quiescent, normal galaxies than theirs. A full 
treatment including the upper-limit sample will be pre- 
sented elsewhere (Takeuchi et al. 2004, in preparation). 
Our final subsample contains 1420 galaxies. 

3.2. Results 

3.2.1. The L(12 /im)-LiR, relation 

For the 12 ^m luminosity, we obtained the regression pa- 
rameters for Lir as follows: 



logiiR = 1.02 + 0.972 logL(12 /tm) , 
r = 0.976 , 
V r = 0.0238 , 



(6) 
(7) 
(8) 



where r is the correlation coefficient and V e is the dis- 
persion in the linear model (see Appendix B). Here the 
above V c gives the 95 % confidence interval 0.3-0.4. The 
data points and the regression line are shown in Fig. [21 
The 95 % confidence limits for the prediction error are 
presented by dotted lines. 

We see a tight linear relation between L(12 fim) and 
Lir, with a correlation coefficient r = 0.976. As seen in 
Sect. |21 the scatter in Fig. [21 is not due to the estimation 
error, but is caused by the intrinsic properties of individ- 
ual galaxies: it is a reflection of the physical variety in the 



Fig. 3. The relation between 25 /im luminosity L(25/im) 
and Lir. 

SEDs of the sample galaxies. We will discuss the origin of 
the scatter in future work (Takeuchi et al. 2004 in prepa- 
ration) . It gives the prediction error of a factor of 4-5 at 
the IR luminosity range [1O 6 L0, 10 n L©]. It is an interest- 
ing result because we know there is a large variety of IR 
SEDs among galaxies, depending on their activities. 

3.2.2. The L(25/im)-L IR relation 

As above, for the 25 /xm luminosity, we obtained the re- 
gression parameters for Lir as 



log Lir = 2.01 + 0.878 log L(25/tm) , 
r = 0.980 , 
V = 0.0203. 



(9) 
(10) 

(11) 



This yields the 95 % confidence interval of 0.3-0.5. The 
data points and the regression lines are shown in Fig. |21 
Again, the 95 % confidence limits are presented by dotted 
lines. The width of the confidence interval corresponds to 
a factor of 4-6. 

Thus, we conclude that both L(12/tm) and L(25/mi) 
provide us with reliable estimates for the total IR lumi- 
nosity Lir, which are valid for several orders of magnitude 
in IR luminosity. 

4. Discussion 

4.1. Applicability and limitation of the linear relations 

In Sect. |2| we obtained fairly tight linear relations between 
MIR luminosities L(12/im) and L(25/mi), and Lir. We 
also found that the scatter in the relations is due to the 
intrinsic properties of the SEDs of galaxies, and we see 
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some galaxies significantly deviating from the 95 % confi- 
dence intervals. Then, a natural question is: for which type 
of galaxy does the relation work well? Among the sample 
galaxies in Tabled we have some galaxies with SEDs in- 
dicative of warm or hot dust (SBS 0335-052, II Zw 40, 
and IRAS F10214+4724), as well as those with SEDs in- 
dicative of cold dust (NGC 1569 and Arp 220). In order 
to examine the applicability and limitation of the rela- 
tions, we revisit the well-known galaxy sample presented 
in Table ^ We represent the luminosity predicted from 
the linear relations [Eqs. 10 and @] by L^ oar . 

We plot the relation between the true integrated L\^ al 
and L^ r /L\^ al in Figs.Han dEl We also show the direct 
estimates from the formula of ISanders fc Mirabell l)l996|) 
using the four IRAS flux densities (filled squares). The 
ratios i}'^ car / £iR tal are presented by open squares with 
error bars that represent the 95 % confidence interval. In 
Fig.0]the prediction is obtained from the 12-/xm relation, 
while in Fig. it is obtained from the 25-^tm relation. 

In Fig. 01 most of the normal galaxies give reason- 
able agreement between Lir and the estimates from the 
linear relation, L\^ clii '. However, the linear relation un- 
derestimates luminosities for three IR luminous galax- 
ies (NGC 2623, UGC 8387, and Arp 220). We also find 
that L^ ear of NGC 1569 is also smaller than the true 
value. In fact, they have strongly extincted, red SEDs (see 
Appendix A), i.e., it is more IR- luminous than expected 
from their MIR luminosities. For the other extreme, the 
linear relation gives acceptable estimates (SBS 0335—052 
and IRAS F10214+4724) within the 95 % confidence level. 
Thus, we conclude that the 12-//m linear relation can be 
applicable for most of the variety of SEDs, except the 
extremely extinguished ones like Arp 220. For such 'red' 
galaxies, it gives a significant underestimation for Ltr. 

In Fig. |SJ in contrast, Arp 220 and other red galax- 
ies are no longer serious outliers. On the other hand, 
SBS 0335—052 significantly deviates upward from the 
true I/tr. Since SBS 335—052 has very hot dust emis- 
sion (Dal e et alJl2001b|) . the linear relation overestimates 
the Lir. Anther two dwarf galaxies, NGC 1569 and 
II Zw 40, are also fairly overestimated because they 
also have warm dust emission. However, the estimate for 
IRAS F10214+4724 is excellent. Hence, the linear relation 
between L(25/xm) and Lir tends to overestimate the Ltr 
for the galaxies with hot dust, but it works well for AGN- 
like SEDs, i.e., SEDs with a hot dust emission as well as 
with a FIR thermal emission. 



10.0 



12/i,m 



4.2. Formula for galaxies at z ~ 0.6 based on 15-fim 
luminosity 

Now we consider the higher-z universe. As mentioned in 
Sect. H our relations will be undoubtedly useful to esti- 
mate Lir for galaxies detected in the very deep Spitzcr 
MIPS data. 

For galaxies at z = 1, the L(12 /j,ul)-Lir linear relation 
itself can be used as an estimator of the total IR luminos- 
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Fig. 4. The performance of the linear estimation formula 
obtained from the L(12/im)-LiR relation. Filled squares 
are the same as those in Fig. ^ The open squares are 
the estimates L} I R lear , obtained from the linear relation 
[Eq. JJjJ], normalized to Lir. The vertical error bars corre- 
spond to the 95 % confidence interval shown in Fig. El The 
linear relation gives reasonable values for normal galaxies. 



ity from the MIPS 24 /im band. What should we do to 
estimate the total IR luminosity for galaxies at redshifts 
between and 1? In this subsection, we try to make a 
useful 'interpolation' formula, which can be used to esti- 
mate the total IR luminosity for galaxies at z ~ 0.5-0.7 
in Spitzer data. 

The practical difficulty is the complexity of the MIR 
SED of galaxies. At these w avelengths, we observe many 
aromatic band features C e. g. . iMaddenl EoOO ) . thus, simple 
linear interpolation might not work well. A more complex 
and continuous interpolation requires some kind of galaxy 
SED model which is no longer free of assumptions, often 
not well-understood. In this work, we stick to the empiri- 
cal relationships directly obtained from observed datasets. 
Thus, based on the ISO deep 15-/xm observations, we try 
to find a relationship between 15 /itm and £ir, since the 
observed wavelength of 24 fj,m corresponds to the emitted 
wavelength of 15 /im at z = 0.6. Although it cannot cover 
the whole range of z = 0-1, it can be applied to a signifi- 
cant fraction of galaxies in this redshift range: taking into 
account the photometric redshift uncertainty, we consider 
galaxies at z ~ 0.5-0.7. If we suppose a flux density limit 
of 200 fiJy, the corresponding luminosit y at these redshift s 
will be L(lbfJ,m) ~ 1O 9 L ~ L* (see iFang et alJll99l . 



Hence, the fraction of galaxies at z 
detected galaxies will be 20-40 %. 



0.5-0.7 among the 
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Fig. 5. The performance of the linear estimation formula 
obtained from the L(25 /im)-LiR relation. Symbols are as 
in Fig. except that they are for 25 /im. 

4.2.1. Estimation formula for the total IR luminosity 
from lb-fim luminosity 

bale et all (|200l4 provided average flux density ratios 
for IRAS and ISO bands as a function of the ratio 
Seo/SioQ. It is well known that these flux density ratios 
depend on the <S6o/<Sioo ratio in gener al, so that the em- 
pirical SEP models wo r k well (e.g.. iDale et al.1 l2001eJ: 
Franceschini et alJ l200lt jTakeuchi et allbOOlat IXu et alJ 



200lHTotani fc Takeuchifl2002HLagache " et al.ll2003|) . For 

our purposes, however, the S12/S15 ratio only weakly de- 
pends on the / S100 ra tio compared to other wavebands, 
because the wavelength difference of these two bands is 
small. We can also derive the formula for 15 /im from the 
L(25/im)-LiR relation via S15/S25, however S15/S25 has 
a stronger and more systematic dependence. Since such 
a systematic dependence will result in a larger dispersion 
in the linear relation and reduce its reliability, we adopt 
S12/S15 for further discussion. 

Then, considering the error of this ratio, we can 
safely use the a verage value over the sample of 
IDale et all l)2001a(l (their Table 1, column 8). We 
found log (Si 2 / Si 5) = 0.112, which corresponds to 
log[ J L(12/im)/L(15^m)] = 0.209. Assuming that the 
slope of the MIR-total IR luminosity relation does not 
change significantly between 12 and 15 /im, we obtain the 
following relation 



logL m = 1.23 + 0.972 log L (15 /urn) 



(12) 



The linear formula between L(15 /zm)-LiR luminosities 
[Eq. (JT^Jl] shows a good agreement with the relation by d i- 
rect fitting of the data proposed bv lCharv fc Elbadll200ll) . 



within the quoted error: 

logL m = (1.05 ±0.174) + 0.998 logL(15/xm). (13) 

4.2.2. Examination of the 15 fim formula by observed 
galaxy sample 

In order to check the validity of Eq. i|12|) . we use 
the quiescent galaxy sample in the Virgo cluster and 
the Coma/Abell 1367 sup ercluster regions ijBoselli et alJ 
|2003|) . iBoselli et alJ (|2p04) have reported a good correla- 
tion between L(15/im) and Lfir for the galaxies in the 
sample. We again constructed a 'good quality' subsample 
with flux densities in all the bands of IRAS and ISO. We 
put a further c onstraint that the d etected flux has qual- 
ity flag 1 [Q of IBoselli etaH (|2003fk column (14) in their 
Table 2] and examined if the flux density suffers contam- 
ination by their close neighbors, and end up with a final 
subsample of 32 galaxies. 

We plot this sample and our empirical formula (with 
95 % confidence interval) in Fig. The formula is rep- 
resented by the solid lines, and the confidence limits are 
shown by dotted lines. Indeed, 31 out of 32 galaxies lie in 
the confidence interval in each panel, i.e., the prediction 
from the formulae successfully work for ~ 95 % of the sam- 
ple. Thus, we conclude that Eq. ^|is a reliable estimator 
of the Lih from 15-/im luminosity with an uncertainty of 
a factor of 4-5, and if the effect of the evolution is small, 
this relation can be used as an estimator of Lju from the 
24 /im luminosity of a galaxy at z ~ 0.6. 

However, we must keep in mind that there is 
clear evidence of a strong evolut ion of galaxies (e.g., 
iTakeuchi etal]l2000L l2001al l2003b|) at < z < 1, and 
we expect a significant brightening of g alaxies up to a 
factor of a few at z = 0. 5 — 0.6 (e.g., ITakeuchi et alJ 
l2001at lLagache et all [20031. Further investigation with 
physically-based models and high-z observations should 
be done in order to examine and/or modify the present 
formulae. 

5. Summary and conclusion 

In this wor k, we first examined four IR luminosity estima - 
tors, L tttr faelou et alJll988fl . L T ir (bale et al.ll2001al) 
£ T tr 2 l|Dale fc Heloull2002|) and L m l|Sanders fc Mirabel 
1996) with the observed SEDs of well-known galaxies. We 
found that £tir, ^tir2, and Lir correct the contribu- 
tion from the wavelengths missed by Lfir, but the lat- 
ter two are better. The estimator Lir provides excellent 
estimates for a very wide variety of galaxy SEDs, from 
SEDs indicative of very hot dust (e.g., SBS 0335-052 and 
IRAS F10214+4724) to very extinguished SEDs and/or 
cold dust emission (e.g., Arp 220). We also note that the 
performance of Ltir2 is almost as good as that of Lir. 

Using Lir, we t hen statistically anal yzed the IRAS 
PSCz galaxy sample l|Saunders et al.l2000jl and found use- 
ful formulae relating the MIR monochromatic luminosi- 
ties [L(12 fim) and L(25/xm)], and Lir. For this purpose 
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Fig. 6. The relation between 15 /im luminosity and L tr for 
the quiescent galaxy sample of iBoselli et alJ (12003^ with 



our formula for estimating Lir. We also show the 95 % 
confidence interval, which is the same as that obtained for 
the L(12/im)-LiR relation. 



we constructed a subsample of 1420 galaxies with all four 
IRAS band (12, 25, 60, and 100 fim) flux densities. We 
found linear relations between Ltr and MIR luminosities, 
L(12/im) and L(25/im). The prediction error with 95 % 
confidence level is a factor of 4-5. Hence, these formulae 
are useful for the estimation of the total IR luminosity 
^totai on iy f rom 12 yjm or 25 /im observations. 

We further tried to make an 'interpolation' formula 
for galaxies in the middle of z — and 1. For this pur- 
pose we construct the formula of the relation between 15- 
/im luminosity and the total IR lu minosity using the flux 
density ratio of lDale et alJ l)2001a(l . The obtained formula 
well reproduced th e observed relation in the sample of 
IBoselli et ai1l|2003|) . We conclude that the 15 -/im formula 
can be used as an estimator of the total IR luminosity 
from 24 /im observations of galaxies at z ~ 0.6. 



Appendix A: SED of our well-known galaxy 
sample 

In Appendix^] we present all the observed SEDs of galax- 
ies we used in examining the performance of the total IR 
luminosity estimators. We show the normal galaxy sample 
with 10 10 L Q < L^ tal < 10 11 Lq in Fig.EU IR-luminous 
sample in Fig. IA.2| and dwarf sample in Fig. lA.3| Among 
the dwarf sample, for SBS 0335-052 and II Zw 40, the 
interpolated points are represented by filled squares (see 
main text). 
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Fig. A.l. The galaxy sample with total IR luminosities, 
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Appendix B: Regression analysis 

We made a regression analysis for the logarithms of L(X). 
It should be noted here that we are interested in estimat- 
ing the total IR luminosity L\^ al from the MIR luminosity. 
Then, in the regression analysis, the uncertainty that we 
need is the so-called prediction error, not the error of the 
regression parameters. We represent the linear regression 
model as 



y = (3 + Pix + e 



where 



iV(0,a 2 ) 



(B.l) 



(B.2) 



Here the symbol l ~' means that the stochastic variable e 
obeys a Gaussian distribution N with a mean and dis- 
persion a 2 . The following estimators bo and b\ are known 
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Fig. A. 3. The dwarf galaxy sample, with L 
1O 1O L0. No observational data are available at FIR- 
submm for II Zw 40 and SBS 0335—052, we extrapolated 
the SEDs by the SED model for a young galaxy presented 
by Takeuchi et al. (2003). The extrapolated data are rep- 
resented by filled squares. 



as the best unbiased estimators 3 for (3q and f3%: 

EiUte -x){vi -v) 



(B.3) 



3 That is, E [6o] = A) and E [6i] = /3i, where E [x] represents 
the expectation value of a stochastic variable x, and the vari- 
ance is the smallest among the estimators. 



and 
bo = V 



(B.4) 



where x and y are the sample mean of Xi and yi, respec- 
tively. The dispersion of bo and b\ shows the statistical 
uncertainty of parameters. However, in a practical appli- 
cation, we need a dispersion of the estimation value y(xo) 
(here hat means that the value is the predicted one and 
not the sample value which would be obtained in a po- 
tential new observation at Xq) for a certain value Xq of 
the independent variable x, in the sense that if we could 
repeat an observation n times, we want an interval within 
which, for example, 95 % of the prediction values y(xo) lie. 
This range is the prediction error, and can be evaluated 
by the formula 



Vr, 



v[y(xo) -y(xo)} 

V[y(x )}+V[y(x )} 

(X ~ xf 



1 + 



(B.5) 



where the symbol V [x] signifies the variance of a stochas- 
tic variable x. The second line of eq. (|B.5|I follows from 
the statistical independence between y(xo) and y(xo). 
Observationally, we should replace a 2 with its unbiased 
estimator, V c [whose unit is a square of dex (an order of 
magnitude in luminosity)], obtained as 



V B = 



(B.6) 



The 95 % confidence interval for the regression line is, 
then, represented by y(xo) ± 2.228-^/V^. For further sta- 
tistical discussions, see e.g., Stuart, Ord, & Arnold (1999). 
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